Maximum Entanglement in Squeezed Boson and Fermion States 
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A class of squeezed boson and fermion states is studied with particular emphasis on the nature of entangle- 
ment. We first investigate the case of bosons, considering two-mode squeezed states. Then we construct the 
fermion version to show that such states are maximum entangled, for both bosons and fermions. To achieve 
these results, we demonstrate some relations involving squeezed boson states. The generalization to the case of 
fermions is made by using Grassmann variables. 
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Boson coherent states have played an important role in the 
production and detection of Bose-Einstein condensate. The 
experimental tools of laser cooling, magnetic and magneto- 
optic traps have advanced tremendously [1, 2]. This has lead 
to a consideration of producing a degenerate Fermi gas as well 
as condensates of rare isotopes, a fact that has been achieved 
experimentally QSHElSl- These results suggest that a care- 
ful study of the Fermi coherent state and Fermi density oper- 
ators is needed. Some of these concepts have been introduced 
for fermion systems long ago (III S S [ISl^ but only recently, 
Cahill and Glauber [If] have discussed notions such as P- 
function, Q-function and Wigner function for fermions; all of 
them are described as a counterpart of the Boson systems and 
are made possible through the use of Grassmann variables. 

In addition to the idea of achieving a Fermi degenerate gas, 
there is a great deal of interest in studying entangled states 
of multipartite systems lfl2l [Tsll . The existence of entangled 
states is directly related to the nature of the quantum mechan- 
ics formalism, based on the structure of Hilbert space and 
the superposition principle. The present interest in entangled 
states is strongly driven by the understanding that this is the 
focal point of studies leading to teleporting of quantum states, 
from one locus to another, which is also the basic ingredi- 
ent of quantum computers lITil [Tsll . The necessity to study 
such entangled states is greatly increased with the suggestion 
that the conditions for teleporting, however, require specific 
states characterized by maximum entanglement. Measure of 
entanglement has been discussed in the literature in differ- 
ent ways in [13, [H [M mil, and the identification and 
construction of maximum entangled states have also been ad- 
dressed inm. 

Consider the state |V')a6 of a bipartite system (A, B). De- 
fine the reduced density operator pa, by 

Pa = Tr6(|^)Qbba('0l), (1) 

where Trft stands for the trace over the variables of the subsys- 
tem B. A measure of the amount of entanglement of the state 
\'ip)ab is given by the von Neumann entropy associated with 
the reduced density operator f^J], Sa = — Tr(pa In pa). The 
question addressed in this paper is concerned with the maxi- 
mum entangled state within the set of all pure states having a 
given (fixed) reduced energy. 



The entropy Sa is a homogeneous function of first degree in 
its dependency on £a, the energy of the system A. Then since 
we wish to maximize 5a [pa] we require (55a [pa] = under 
the constraints 

£a = {Ha)=TY{paHa), Trpa = l, (2) 

where Ha is the energy operator of the A-component. Follow- 
ing standard procedures, we derive then a constraint equation 
for Pa; that is Aq — 1 + XiHa — In pa = 0, where Aq and 
Ai are Lagrange multipliers. Solving this equation, we get a 
Gibbs-like density operator, that is 

Pa = -^exp(AiiJa), (3) 

where Z — exp(l — Aq). Multiplying the constraint equa- 
tion by Pa, taking the trace and using Eqs. (|2|l and ((Sj, we get 

\nZ = Alfa + Sa- 

For the sake of convenience, let us write Ai = — t, then 
we have —r^^lnZ = £a — T~^Sa- The function F{t) = 
— r^^ In Z describes the Legendre transform of Sa since we 
assume that t = dSa/d£a- Here r is an intensive parame- 
ter describing the fact that the energy average £a = {Ha), 
given by Eq. Q, is constant in the state described by pa- 
Therefore, a state |V-')ab with reduced energy £a is a maxi- 
mum entangled state when the corresponding reduced density 
matrix, defined in Eq. ([T]i, is written as a canonical Gibbs- 
ensemble distribution, explicitly given by Eq. (O, such that 
Z = [1 — exp(— r)]~^. Using this approach, examples of 
maximum entanglement states were explicitly constructed in 
ll23ll . considering as a guide the thermofield dynamics (TFD) 
formalism. 

In the present work we extend that analysis to show that 
two-mode squeezed states (TMSS) are maximum entangled, 
within a class of states with fixed energy, for both fermions 
and bosons. We first investigate the case of bosons and, then, 
we construct the fermion version, extending a preliminary 
study by Chaturvedi et al. 1 10,,]. For the case of fermions, the 
situation is more intricate, demanding the notion of coherent 
fermion state and density operator, which is achieved by using 
Grassmann variables. In any case we explore the similarities 
among these squeezed states with those found in the TFD for- 
malism iH HI [M]. 
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In order to analyze entanglement, we derive some gen- 
eral pro pert ies of the usual Caves-Schumaker (CS) state for 
bosons 1I27I1 . Consider a two-boson system specified by the 
operators a and b obeying the algebra [a, a^ = [b, b^ = 1, 
[a,b] — 0, with unitary displacement operators -Da(C) = 
exp[^a^ ~ a^*] and Db{rj) — exp[776^ — brj*], and define the 
two-mode squeezing operator 



Sabil) = exp[7(at5t - ab)], 



(4) 



with 7 being a real non-negative number for simplicity. Using 
standard TFD formulas [ 24. ,25il . we get the useful relations 

0(7) Sab(.l)aSl^,{j) =u{-i)a~v{-i)b\ 
b{i) = Sab{l)bSlf,{-i) ^ u{-i)b- v{-i)a\ 

and the corresponding relations for 0^(7) and ^'^(7), where 
14(7) — cosh7 and 17(7) — sinh7. 

First, consider the TMSS defined by \"f)ab — <S'a6(7)|0)Q6, 
where |0)q6 — |0)q 0)b = |0)a|0)h is the two-mode vac- 
uum such that a\Q)a — fe|0)b = 0. For this squeezed 
state, we have a(7)|7)„h — 6(7)|7)ab — 0. Another 
important result is that Sab{l) is a canonical transforma- 
tion, that is, [a(7),a^(7)] = [^(7)1^^7)] = 1 ™d 
[0(7), 6(7)] = 0. Now, using the operator identity exp[7(^ + 
B)] = exp[(tanh7)i3] exp[(lncosh7)C] exp[(tanh7)74], 
with A = -ab, B = atfct and C = [A, B] = -a^a - bb\ the 
TMSS can be written as 



\l)ab = 



1 



cosh 7 



exp[(tanh7)a'l'6'l'] |0)afc. 



(5) 



Changing the parametrization by taking cosh 7 — [1 — 
exp(— t)]~^/^, so that tanh7 = exp(— t/2), and defining 
Z{t) = [1 — exp(— r)]~^ = Tr exp[— ra^a], we find 



|7>q6 



1 



Therefore, the TMSS \^)ab can be written as 



n)a\n)b- 



\l)ab = \fjjj) X! W)a\n)b, 
n=0 



(6) 



(7) 



where fair) = exp{-'Ta^ a)/Z{T). 

The CS states are introduced by the application of 
Sabil) to a two-mode coherent state, that is |^,?y,7)) = 
Sabil)Da{ODbiv)\0)ab- We can show that |C,?y,7)) = 
le" r?,7), with |e,^,7) = Da{ODb{fl)Sabhmab, where 



Bb{i) being the matrix form associated to Sab{l) given by 



(8) 



Bb{i) = 



"(7) -«(7) 
-w(7) m(7) 



(9) 



For ^ = r/ = 0, the CS state reduces to the TMSS \^)ab = 
•Safe (7)10), which has the same structure as the thermal vac- 
uum state used to introduce TFD ll24i I25I1 . We explore this 



result to show that 77, 7)) can be used to define a Gibbs-like 
density. 

Following the prescription discussed before, taking pab — 
ICj "Hi 7) (7j CI' we have to perform the calculation of the re- 
duced density matrix, say pa — TibPab- Using the notation, 
|r)b = n!^^/^(&^)'"|0)f, (similarly for mode a), we can write 
the matrix elements {s\pa\t) = Y.r '^'\^'^\b{APab\r)b\t) a as 

{s\pa\t) - a{s\Da{0^/Ur)\n)ab{r\D,{j^)\n), 



X b{fn\Dbijl)^\r)ba{m\^/Uj)Da{i)^\t)a. 

Changing the order of the matrix elements in the b mode, 
and using the completeness relation, we obtain {s\pa\t) — 

{s\Da{Ofa{T)Da{S,y\t). ThuS, WC get 

Pa = DaiOfaiT)Di{0 = (t) exp(-™t (^)a(^) ) , 

(10) 

where a(^) = Da{£,)aDl^{£_) is the displaced annihilation op- 
erator; Pa is, therefore, a Gibbs-like density. In particular, for 
^ = 0, we find pa — fa{T) showing that the TMSS |7)q6 also 
generates a Gibbs-like density. 

Using the displaced Fock's basis, {Da{£,)\n)a}, we show 
that the reduced von Neumann entropy for a CS state is 



Sir) 



- 1 



ln(l 



(11) 



Thus, all CS states, with the same (fixed) squeezing param- 
eter, have the same amount of entanglement independent of 
the displacement parameters. Among them, the one having 
the smallest energy is the TMSS (^ — 0); its reduced energy 
{£a = Tr{paa''a)) is given by £{t) — (e"^ — 1)~^. Since 
both actions of displacing and squeezing the vacuum lead to 
states with greater energy, the TMSS is the maximum entan- 
gled state when the energy is fixed. 

Let us now analyze comparatively the amount of en- 
tanglement of the TMSS. Consider the state l^*)^^^ = 

'^n=o 1*^)0 l"')6 which has reduced energy and entropy 
given by £'{N) = {N - l)/2 and S'{N) = InN. This state 
has the greatest amount of entanglement among all pure states 
belonging to the finite (iV^) dimensional subspace spanned by 
{|0). |0), , |0)„ ,...,\N-1)JN- 1) J,coii-esponding 
to equal occupation probability. Naturally, as iV ^ 00, both 
energy and amount of entanglement of goes to 00. 

Now, take another parametrization of the TMSS by writ- 
ing r = ln(x + 1) — ln(x — 1); the limit situations of 
zero and infinite squeezing correspond to x = 1 (7 = 
0, T — 00) and X = 00 (7 = cx), t ~ 0), respec- 
tively. The reduced energy and the amount of entangle- 
ment are then written as £{x) = (x ~ l)/2 and S{x) — 
iix + 1) ln(x + 1) - (x - 1) ln(x - 1) - 2 In 2] /2. We find 
that both £{x) and S{x) go to 00 as x ^ 00, with S{x) ~ 
Inx in this limit. In Fig. 1, we plot the difference between 
S{x) and Inx, showing explicitly that the TMSS has an 
amount of entanglement greater than that of the state 
with the same energy, for all N > 2. 
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FIG. 1 : Plot of AS = 5(x) — In x as a function of x- 



We now consider a similar analysis for fermions. A co- 
herent fermion state (j7i (st liH] can be defined by introduc- 
ing a displacement-like operator D{a), where a is a Grass- 
mann variable, trying to reproduce formally the basic re- 
sults of the boson case. This is accomplished in the follow- 
ing way. Consider a and f3 two Grassmann numbers, then 
{a, /3} = a/3+/3a = 0, and moreover, {a, a} — {a, = 0, 
where the fermion operator a and satisfy the anticom- 
mutation relation {a,at} = 1. Notice that we maintain the 
same notation for the creation and annihilation operators as 
that for bosons. The complex conjugation is an antilinear 
mapping * : a* such that, for a general expression 
involving Grassmann numbers and the operators a and a^, 
we have {aia* + c(3il3j)* = a-jO* + c*l3*j3*; c e C and 



aka\a*fii. The Grassmann 



variable a is considered independent of a* . 

The fermion displacement operator is defined by 

Da{a) = exp(a^a — a*a)^ 



(12) 



such that Da{a)aD'l{a) = a — a, with the coherent state 
given by \a) = Da{a)\Q)a and a\a) = a\a). We can then 
prove that the dual coherent state is given by (a| = {Q\D\{a), 
with -Dj(a) = D'^^{a). As a consequence 



exp(a*/3 — 2'-'^*'^ 



and {a\D\{a) = {a\a* . In terms of the number basis, the 
state la) is written as 



-a' a/2 ' 



(13) 



and, then, we have {n\a) = exp(— a*a/2)(Q!)". 

The integration is defined, as usual, by j da — and 
/ daa = 1. Note that, in particular, we have J da*a* — 
1 (resulting in (da)* = —da*), J da*daaa* — 1 and 
/ da* da \a) {a\ = 1. 

Cahill and Glauber ifTlll introduced the following coherent 
fermion state representation for a density operator, i.e. p = 
/ d^aP(a) I — a){a\, where P{a) is the corresponding P- 
function. Notice that the density operator 



Pa 



a){a\ 



(14) 



possesses the expected properties to be taken as representing 
the coherent state |a). First, it is normalized, Tipa = 1. This 



property can be proved from the matrix representation of pa, 
by calculating {m\pa\n) — exp(aa*)(— Q;)'"(a*)", giving 
rise to 



1 — a*a a* 

—a a* a 



Secondly, TT{paa^a) — a*a, which is similar to the boson 
case. Observe that pa, although not being hermitian, is intro- 
duced in such way that p^ = p fill] . 

Using the properties described above, we can prove that the 
displaced fermion number state is given by 

Da{a)\n)^{a^ -a*r\a), 

with n = 0, 1. Another property useful for calculations but 
reflecting also the nature of pa is given by 



(m|p„|n) = (-ir("+i)(a|n)(m|a), 



(15) 



where {m\pa\n) = (m| — a){a\n). Observe that for m ^ n 
we have {n\ ~a){a\n) = {a\n){n\a). The usefulness of this 
result can be seen in the proof that Ti pa = 1. 

Let us now consider a two-fermion system, specified by the 
operators a and b satisfying the algebra {a, a^} = {b, b^} — 1, 
with all the other anticommutation relations being zero. A 
fermionic two-mode squeezed vacuum state is defined by 
\l)ab = Sab{l)\0)ab, whcre 7 is still a real number and 
Sabil) — exp[7(a^5t — ab)]. Some useful formulas can be 
derived using 5*06(7), that is, 

0(7) = Sab{l)aSl^{-f) ^ u{j)a- v{j)b\ 
6(7) = Sab{l)bSl{j)^u{j)b + v{j)a\ 

where now ^(7) ~ 003(7) and 11(7) — sin(7). Thus, for the 
two-mode squeezed vacuum state \j)ab we have a{j)\'y)ab = 
b{l)h)ab = 0, since a|0), = b\0)b = 0. 

The squeezing operator Sabil) is a canonical trans- 
formation, in the sense that, as in the case of bosons, 
{a(7),a(7)t} = {5(7),6(7)t} = 1 and {ai-/),b{j)} = 0. 
The matrix form Bpij) associated to Sabil) is 



"(7) 



(16) 



The vector \'^)ab can be cast in a TFD state. To see that let 
us write 

b)ab = [l-^{ba-a^b^) + ^{ba-a^b^)^ + ...]\0)ab. (17) 

Using the relations {ba - atfet)2n|o)^^ = (-l)"|0)ah and 
introducing the reparametrization 11(7) = cos 7 = (1 + 
e^"^)^-^/^, v{j) — sin7 = (1 + e^)^^/'^, we obtain 



l7>afc 



1 



(l + e--/V6t)|0)afc. 



\/T+7 

Defining Z{t) = 1 + e"^, Eq. ^ reads 
1 



(18) 



h)ab = 



e-'''^H\0)a\0)b + \l)a\l)b), (19) 
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where N = a^a, the fermion number operator for the mode a, 
is such that N\n)a — "|?^)a• Therefore, we obtain 

1 

n=0 

with /a(T) = Z^^{t) exp(— Tfl'l'a). 

With these resuhs, we can prove the following statement. 
Given the two fermion displacement operators, Da{a) = 
exp(a^a — a*a) and Dh{(i) — exp(&^/3 — where a and 
/3 are Grassmann numbers, then 

Sab{l)Da{a)Di{(3) = Da{a)Di{P)Sab{l) (20) 

where 

(ji)=M7)(;.), (21) 

Thus, the fermion version of the CS state, defined by 

\a,P;^)) = Sab{l)Da{a)DbmQ)ab, 

is related to the state |a,/3;7) = £'Q(a)£'b(/3)5'Q6(7)|0)afc 
by the transformation given in Eqs. ( |20l l and ( 1211 1. As in the 
bosonic case, when a = /3 = we have the two-fermion 
squeezed vacuum state \-f)ab = <S'ab(7)|0)a6. 

Now we turn our attention to the nature of the entanglement 
in squeezed fermion states. Considering the states |q;,/3;7) 
and inspired by the definition of the density operator given in 
Eq. (IT4b . we introduce the following density matrix 

Pab = I -a,-/3,7)(7,^,a|. 



Performing the trace in the mode h and using the properties 
derived before, we can prove that pa = Da{a)fa{T)Da{ay, 
similar to the boson case. Thus we find that the state |a, /?, 7) 
has reduced density operator in the form of a Gibbs-like den- 
sity. The reduced entropy is thus maximal. However, in 
the fermionic case, the CS states are not in general physical 
states fTT] since they involve Grassmann variables. The two- 
fermion squeezed vacuum state \ j)ab is nevertheless physical 
and maximally entangled. 

It is worth mention that, in the case of fermions, there is an- 
other class of physical states having maximum entanglement 
for a given value of the reduced energy. In fact, one can show 
that the state 

bYab = {\0)a\l)b + e--/^|l)„|0),) / (22) 

has reduced density operator p'^ identical to the reduced den- 
sity operator pa associated with the state \'y)ab', therefore, 
these states have identical reduced energy and entropy. 

Summarizing, in this paper we have analyzed a class of 
two-mode squeezed boson and fermion states, looking for 
explicit realization of maximum entangled states with fixed 
energy. We investigate the case of bosons, and then, con- 
struct the fermion version, to show that such states, in both 
cases, are maximum entangled. For achieving these results 
we have demonstrated some relations involving the squeezed 
boson states, which are then extended to the case of fermions. 
The calculations for fermions are performed with a general- 
ization of the density fermion operator introduced by Cahill 
and Glauber (u^. 

This work was partially supported by CNPq (Brazil) and 
NSERC (Canada). 
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